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Abstract 

Recently, VT symmetry of many single-particle non-Hermitian Hamiltonians has 
been conjectured sufficient for keeping their spectrum real. We show that and how 
the similar concept of a "weakened Hermiticity" can be extended to some exactly 
solvable two- and three-particle models. 
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1 Introduction 



Insight into bound states in quantum mechanics is facilitated by solvable models. 
They clarify the structure of the single-particle wave functions in the context of 
supersymmetry Lie algebras || and Sturm-Liouville oscillation theorems ||. 
This approach can immediately be extended to the systems of more particles where 
a very exceptional role is played by the Calogero's exactly solvable Hamiltonian 
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i=l Udj i- i<j=2 
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It describes A particles on a line with 2m — H — 1 and its Lie algebraic treatment 
proves enormously productive [H. The analysis and interpretation of some of its 
properties is also facilitated in the non-singular limit g — > where the interaction 
degenerates to the mere harmonic-oscillator attraction. 

Recently, the so called VT symmetric quantum mechanics of Bender et al |J 
offered a new picture of some single-particle models. For example, the most common 
spectrum of the harmonic oscillator in D dimensions was re-interpreted as a special 
case of a non-equidistant real spectrum pertaining to a slightly more general non- 
Hermitian model ]/]]. 

The picture is based on a complexification of coordinates which breaks the Her- 
miticity of the Hamiltonian but does not destroy the reality of the energies. In the 
present paper we intend to show that such a complexification method can be gen- 
eralized and applied to some many-particle Hamiltonians. In a constructive way we 
are going to demonstrate the existence of non-Hermitian modifications of and 
if( 3 ) which consequently preserve the reality of the spectrum. 

2 Acceptable solutions 

At any integer A > 1 the introduction of the centre of mass 
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enables us to eliminate the bulk motion of the Calogero system. Submerging it for 
this purpose in an external auxiliary well (R^) = Aoj 2 [R^] 2 /8 we get a slightly 
simplified version of the Calogero Hamiltonian, 



i=l 



d 2 A 2 2 

+ — LU 2 xf 



dx 



A 

+ E 

i<j=2 



\Xi Xjf 



This leads to the centre-of-mass equation 



+ -u 2 R 2 -u \ -F 



dR 2 8 



V(R) = 



with the well known solutions pi. Each element of its spectrum 



N 



2N + 1, 



N = 0,1,... 



adds a constant to all the internal energies E. 

In the first nontrivial model with the mere two interacting particles let us 
ignore the centre-of-mass equation (|l|) as trivial and skip g = (harmonic oscillator) . 
Our Schrodinger bound state problem is then described by the singular ordinary 
differential equation in the relative coordinate X = (x-y — x 2 )/\/2, 



d 2 
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d X 2 + r 2X2+ 2± 2 - E 



= o. 



(2) 



The symbol E denotes the energy in the centre-of-mass system. The singularity 
resembles the current centrifugal term with a non-integer real parameter 
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(3) 



Such a singularity is, mathematically speaking, too strong. One has to reduce the 
domain to a half-axis and fix the ordering of the particles (say, x\ > x 2 , i.e., X > 0; 
sometimes, this is interpreted as a choice of the Bolzmann statistics). 

In a mathematically more rigorous setting one should even demand that the 
repulsion is not weak, g > 3/2. Otherwise, the Hamiltonian admits many self-adjoint 
extensions, each of which may lead to a different spectrum 0, [H^ . 

In the Calogero's paper HJ this problem has been addressed and resolved in a 
pragmatic spirit. Even when the irregular solution ip^~'(X) ~ |^| - ^ \X\ 1 
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becomes normalizable (which certainly happens for all t < 1/2 and/or g < 3/2) 
we eliminate it as "physically inacceptable" via an ad hoc condition. The details of 



this argument may be found in the recent comment |L0[ where the eligible ad hoc 
conditions were listed as depending on the range of g, 

lim M [X- 1 ^ Hermttia ^{X)] = 0, g E (-1/2,0), 

lim M ^(HermUian)^ = Qj g g (0,3/2), (4) 

lim x ^ ^(Hermitian)^ = Qj g £ [3/2, oo). 

The constraint g/2 = l{t + 1) > —1/4 in the first line is unavoidable and protects 
the system against a collapse |IT| , while only the third line is fully equivalent to the 
conventional requirement of normalizability |~2] . 



In what follows we shall use the same philosophy. The validity of the appropriate 
conditions of the type (f|) will be postulated as a conventional regularization of the 
singularities. In the spirit of ref. 0] this will define the "acceptable" solutions and/or 
make them unique. 

In the well known A = 2 case this implies the termination of the confluent 
hypergeometric series to the Laguerre polynomials, 

^f^^^-X^exp^-^X^L^Q^X 2 ), n = 0,l,... . (5) 

Although the two intervals X > and X < are impenetrably separated, one can 
deal with the presence of the singularity by another physically motivated postulate 

^(bosonic) f—X) = +ll)( b ° SOniC \X) , ^(/ e "™°™ c )(_X) = — ip(f ermionic ) (X) (6) 

mimicking the Bose or Fermi statistics. Such a freedom is rendered possible by the 
full separation of the domains X > and X < 0. By construction, the bosonic and 
fermionic wave functions vanish at the matching point X = 0. In what follows, a 
different type of the ad hoc symmetry will be employed and postulated. The resulting 
form of the modified statistics will be closely related to the so called VT symmetry 
in the single-particle quantum mechanics 0. 
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3 VT— symmetric quantum mechanics 



3.1 A = 1 and the complexification 

No spikes exist in the trivial Hamiltonian It is equivalent to the harmonic 

oscillator but offers still a fairly nontrivial methodical lesson. In their pioneering 
letter ]13| Bender and Boettcher emphasized that the complex shift of coordinates 
R = R(r) = r — ie, r G (— oo, oo) in the harmonic oscillator Schrodinger equation 
([I]) does not change its validity and preserves the normalizability of its wave func- 
tions ^(-R). The spectrum remains unchanged even when we admit a symmetric 
r— dependence in e = e(r 2 ). 

Empirically, the similar coexistence of the real spectrum with the non-Hermitian 
Hamiltonian has been detected for the various other single-particle potentials |IJ] . 
The phenomenon has attracted a lot of attention in the literature |TJ| . Its appealing 
interpretation has been conjectured in terms of the symbols V (which denotes parity, 
VRV = —R) and T (this is "time reversal" or complex conjugation, TiT = —i). 
For many Hamiltonians which commute with the product VT, people have observed 
the reality of the spectra and attributed it expressis verbis to the VT symmetry 



EPS li 



In the present paper we intend to extend this language to cover also some systems 
of more particles. 



3.2 A = 2 and the regularization 

The complex shift of coordinates did not change the solutions of A = 1 equation ([!]) 

in the asymptotic region. The same is true for wave functions of the two Calogero 

particles. In the A = 2 equation (0) the formula 

1 1 2e 2 2ixe . 2 

X2 ~ x 2 + e 2 ~ (x 2 + e 2 ) 2 + {x 2 + e 2 ) 2 ' ^ —x-ze{x ) 

indicates how the deformation of the integration contour regularizes the singularity. 
An immediate practical compensation of the loss of the Hermiticity of is found 
in an improvement of its regularity. This eliminates the rude constraint (|j) and 
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opens space for new solutions. At the same time, we should not get too many of 
them [18] and so we demand the PT symmetry of the contour, 



VT ■ X{x) - VT = X{-x). (7) 

It is possible to search for the even-parity-like solutions of the complex, regularized 
radial equation (§) in a way proposed originally by Buslaev and Grecchi and 



dictated by the angular- momentum interpretation of the parameter £ [19[]. In paper 
we succeeded in complementing the known PT-symmetrized solutions (^) of eq. 
(@) by the second hierarchy, 



^™\ X ) ~ X^exp (-^X 2 ) L-J-^ QcX 2 



71 = 0,1 



At almost all values of g the pair of solutions and (Q) re-connects the subdomains 
x > and x < 0. The equidistance of the old Hermitian spectrum is manifestly 
broken by the new even- like energies. Figure 1 illustrates the result. It displays both 
the energies 

E (Hermttian) = E (+) = I ^ (4 n + 2£ + 3) , E^ ew) = = - U (An - 2i + 1) 

(with n — 0, 1, . . .) as functions of the coupling g (or rather £ = £(g)) at a fixed value 
of the spring constant uj = 2 ■ \J2/A. 

We may summarize that at A = 2 the imaginary shifts of the single-particle 
coordinates x\ and X2 complexify the Jacobi coordinates R and X. This has com- 
paratively trivial consequences for the regular centre of mass problem (p. The 
change proves much more influential in the singular radial equation (fj). In the new 
language we were able to postulate the two types of behaviour near X = 0. Vice 
versa, the current return to the Hermitian constraint (|j) acquires now an unconven- 
tional understanding of the reduction of the spectrum, caused by the pragmatic ad 
hoc elimination of all the redundant states (|). 



4 Generalized statistics 
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4.1 A = 2 and the idea 

The consistent complexification of our two-body singular equation (H) requires an 
introduction of an upward cut in the complex plane of X. This reflects the branching 
role of the singularity and implies that we have to choose e(0) = —ImX(0) > 0. 
Still, having the sub-intervals x > and x < inter-connected, the freedom (|6]) of 
the choice of statistics seems lost. In fact, it is not. This is to be shown below. 

The analysis is facilitated by the symmetry (|7|). We can return to Figure 1 and 
notice that the levels cross at g = 2 k 2 — 1/2 for k = 0, 1, 2, . . .. Everywhere off these 
exceptional points we can start from the g = or £ = bound states with the well 
defined values of parity = ±1. We propose to use a smooth continuation in g. This 
transfers the label ±1 to almost all the energies. We re- introduce the complexified 
bosonic and fermionic symmetry in one of the most natural ways. 

In a less intuitive setting we have to imagine that the parity is not conserved. At 
A = 2 we are able to replace this concept by the equivalent permutation symmetry. 
Our knowledge of the explicit wave functions (^) and (|8|) enables us to speak about 
the bosons and fermions defined by the following rule 

^(bosonic) (—X) = ( — l)~ i ^ bosonic \X), 
^(fermionic) (—X) = ( — \y+ l ermionic ) (X) . 

This is a natural generalization of the I = ad hoc conditions (|6|) prescribing the 
behaviour near the singularity. Mathematically, it enables us to discretize the bound 
state spectrum in a way employed also in refs. 0, [16|, [20 1. 



Whenever we choose an integer i the role of the parity becomes partially re- 
established. The new definition of the statistics (|9|) preserves many features of 
its Hermitian predecessor (^) also in the limit linx^oo e(x 2 ) = 0. The number of 
nodal zeros can be shared by our new bosons and fermions. Their spectra are 
mutually shifted. Geometrically they are formed by the opposite branches of the 
energy parabolas in Figure 1. In terms of an abbreviation a = I + 1/2 the tran- 
sition from fermions to bosons is a mere change of the sign at this square root, 



a = Jl/A + g/2 - -a = -Jl/A + g/2. 
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4.2 A = 3 and a toy model 



The solution of the Calogero three-body problem proves significantly facilitated after 
one omits two of its spikes. For the purely Hermitian toy Hamiltonian 

9 



i=l 



d 2 3 2 2 



<9x 



+ 



(xi - x 2 ) 2 



this has been noted in the first half of Section 3 in ref. |4j]. After a constant imaginary 
shift of the variables xi, x 2 and x% in the three-body Hamiltonian the centre of 
mass R— (xi + x 2 + xs)/y/3 becomes complex. Its elimination is still trivial. The 
other two Jacobi coordinates 

Xi — X 2 xr X\ J rX 2 — 2 £3 



X 



Y 



V2 ' Vq 

enter the partial differential Schrodinger equation 



d 2 



d 2 



■u 2 {X 2 + Y 2 



E 



$(X,Y) = 0. 



dX 2 BY 2 8" ' ' ' ' T 
At A = 3 it replaces the A = 2 radial equation (Q). We employ its separability in the 
hyperspherical coordinates p and <fi. Under the guidance of the Calogeros's paper 
0] this gives X = p sm<f) and Y = p cos</> while 3>(X,Y) = tp(p)f(<j)). The minor 
relevance of the centre-of-mass motion extends naturally to the p— dependence. With 
p E (0, oo) it is described by the wave functions 

3" 



V»n,fc(p) = P m exp 
with the energies 

E = E n h ~- 



'-u[2n+l + (3(k)l 



-LU p 



n, k = 0, 1, . . . 



n, k = 0, 1 



Both depend on g in a way mediated by the new quantity (3{k) > 0. It is defined as 
the square root of the eigenvalue of the "innermost" hyperangular equation 

/ d 2 
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(10) 



\ d<p 2 2 sin 2 y j 

Our key idea is to complexify solely the third coordinate 0. We shall use the recipe 
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inspired by its A = 2 predecessor. The "Bose" or "Fermi" statistics of our toy 
solutions will be determined by their regularized behaviour @ with X(x) replaced by 
0(0- m effect the complex double well problem (|10D is appropriately constrained in 
a way consistent with the permutations of x% and x%, i.e., with the VT-like symmetry 

and with the ordinary real parity conservation reflecting the unconstrained variability 

of x 3 , 

Xk {ix - 0) = (-l) k Xk (<f>) . (11) 
This establishes the closest parallels between A = 2 and A — 3. 

5 Spectra 

5.1 Explicit solutions of the toy model 

Differential equation ([IT]) possesses a general solution which is a hypergeometric 
series of the Gauss type, 

X (±) (0) = (sin 0) 1/2±a 2 F l {u {± \ v (±) ;l± a; sin 2 0), a = -Jl + 2g>0 

with 2u^ = 1/2 — /3 ±a and 2?/ ± -' = 1/2 + /3 ±a and with an equivalent alternative 
form 

x (±)(0) = cos (sin 0) 1 / 2±a 2 F 1 [u {±) + l/2,u (±) + 1/2; 1 ± a;sin 2 0). 

Only the regular, (+)_ superscripted states were acceptable in ref. Q. For us, 
the "irregular" x^K'P) represent bosonic states. These functions are smooth and 
bounded due to our complex regularization of eq. fllPp. We can derive their estimate 
^(±)(0) ~ £ 1 / 2±a [l -f 0(e 2 )] for the non- vanishing e ~ e(0) > and in the closest 
vicinity of the singularities £ = and £ = ±7r, i.e., for £ G (— e, e), £ G (— 7r, — 7r + e) 
or £ G (7r — £, 7r). This implies that both the signs can equally well appear at the 
parameter a. 
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On every boundary of convergence sin 2 = 1 of our power-series solutions we 
arrive at the same necessity of termination as before. In the light of the "real" 



symmetry (|TT|) this occurs if an only if our free parameter belongs to the sequence 



P(k) = (3 {± \k) = k±a + 1/2, 



fc = 0,l, 



In a way paralleling the Calogero's construction this leads to the Gegenbauer poly- 
nomials in 



Xk 



(fermions /bosons) 



(<P) = (sin 4>) 1/2±a ±Q (cos 



n =- A /l + 20>O. 



An analogy of this key formula with its A = 2 predecessor becomes clearer in the 
older notation with 1/2 + a = £ + 1 and 1/2 — a = —£. 

The set of the toy eigenvalues (3 2 decays in the two subsets. They differ just by 
the sign attached to the square-root parameter a. The energy spectrum emerges in 
the compact form 



E 



W = J -u(4n + 2k±2a + 3), a = -Jl + 2g > 0, n,k = 0, 1, 



In the convenient scale JjW — 1, it is sampled in Figure 2. 



5.2 The genuine three-body Calogero-type system 

The resemblance between equations (fj) and (|T0|) is not too surprising since in our 
toy model the third particle is only bound by the purely harmonic forces. Let us 
now recall the full-fledged Calogero Hamiltonian 



d 2 3 2 2 



+ 
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+ 







+ 



(xi-x 2 ) 2 (S2-X3) 2 (a; 3 -a;i) 2 



dxi 2 8 

and recollect that the related partial differential Schrodinger equation 

r) 2 r) 2 S 1 
OX 2 dY 2 8 V ; 2 y 



$(X,F) = 



9 



degenerates to the ordinary differential "innermost" equation 

M / fc (0) = (3 2 (k)M f k {<f>) 

containing the hyperspherical momentum operator 

d 2 1 / o 3 <7i q 2 



1 2 Vsin 2 sin 2 (0 + §7r) sin 2 (0-§7r) 



For the equal couplings Qj = g we have the complex differential equation 

1 + TT^hl) XkW = P 2 (k) Xfc (0). (12) 



V d(j) 2 2 sin 2 3 . / 

It is a straightforward six-well modification of the above toy double-well problem 
with £ G (— 7T, 7r). In its analysis let us start from the Hermitian limit e = with the 
clear geometrical meaning of the permutation symmetry in the real X — Y plane. 
It is divided into six subdomains separated by the impenetrable two-body barriers. 
Thus, the x\ <-» £ 2 interchange Vn-2) happens precisely along the Y — axis. Similarly, 
the other two permutations x 2 <-> £3 and £3 <-> Xi take place along the respective 
lines y = iX/v/3 [§. 

As soon as we introduce e ^ the action of the complexified permutations has 
to be subject to the triple rule 

PT (1 _ 2) -0(O-PT (1 _ 2) = 0(-O, 
PT (2 _3) ■ 0(0 ■ PT (2 _ 3 ) = 0(|vr - 0, ( 13 ) 
PT (3 _i) • 0(0 • 7T( 3 -i) = 0(-|tt - 

which guarantees the commutativity of the Hamiltonian with the vectorial operator 
VT. In analogy with equation (|9|) we postulate that our solutions satisfy the six ad 
hoc boundary conditions 

/( 6osonic )(-^-) = (-l)^/ (bosonic) (0j)> 

0, = 0-|vr, j = -2, -1,0, 1,2, 3. 

This extends the VT — symmetric quantum mechanics of ref. to our complexified 
system of three particles. 
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5.3 Solutions 

The parameter I = £(g) is the same function (§) of g as above. Incidentally, it 
vanishes in the harmonic-oscillator limit g — > 0. Explicit solutions of our Calogero- 
inspired model can be constructed with an ample use of its toy predecessor. In 
particular, the toy hypergeometric solutions have only to be subject to the modified 
boundary conditions. This parallels the Hermitian situation. 

Firstly, due to the overall symmetry (]TJ]) it is still natural to match the wave 
functions in the middle of the separate wells. Formally this relies on the fact that 
all of the infinite hypergeometric series reach their limits (i.e., radii) of convergence 
precisely at/along the three real-middle-of-the-well lines Y — (= X— axis) and 
Y = ±y/3X. For this reason they have to degenerate to the Gegenbauer polynomials 
as before. 

Secondly, the exact values of the energies of the triple-spike complex Calogero 



model follow immediately from a mere comparison of equations (|1QD and fll2|) . The 
appropriate re-scaling of the "intermediate" eigenvalue parameter (3 gives the result- 
ing final form of the spectrum 



Efl = {An + 6k±6a + 5), a = ^-Jl + 2g>0, n, k = 0, 1, . . . . 
V o 2 

It is illustrated in Figure 3. The pertaining complex angular wave functions 

r( f ermionic) / ± \ ( f ermions) / ±\ (■(bosonic) / ,\ (bosons) / nl \ 

fk w) = Xk V>4>), fl '(4>) = Xk (30) , 

and their oscillator-like radial counterparts 

VS(P) = p 3fe±3a+3/2 exp ("\/|-p 2 j ^ ±3a+3/2 

are easily determined using the elementary insertions in the obvious manner resem- 
bling the classical Calogero's analysis. 



'3 _. . _ 1 




6 Summary 

In a way guided by the available one-body experience we complexified the two- 
and three-body Calogero model, having consequently employed the advantage of its 



11 



separability. This enabled us to suppress the majority of complications (e.g., an 
interplay between several complex variables) which would necessarily arise in any 
more complicated many-body model. 

Our conclusion of upmost importance concerns the new role of the operator VT . 
We were guided by its A = 2 action upon X ~ x\ — x 2 as given by equation ([7p. This 
equation was interpreted as a complexification of the usual permutation of particles 
X\ <-> X2- Hence, we also assigned a new meaning to the whole concept of the VT 
symmetry for more particles. Its re-examination in the three-body context enabled 



us to find its appropriate permutation-symmetry generalization (13). 

During our application of the complex permutation symmetry to the Calogero 
model we arrived at several satisfactory analogies between its separate A = 1, A = 2, 
and A = 3 special cases. Firstly, as long as we studied the complex coordinates only 
on the "innermost" level, we were able to leave all the previous stages (with any 
complex Jacobi coordinates) open to virtually arbitrary analytic continuation. 

Secondly, it was comparatively easy to specify the correct physical states in the 
limit of the standard Calogero system. This backward relationship to Hermitian 
predecessors is quite instructive. Among our bosonic as well as fermionic complex 
solutions some of them "survive" in the limit e — > 0. Still, as long as the domains 
themselves become halved, the full set of solutions is overcomplete. For this reason 
the VT symmetric bosons were declared to be redundant. They must be eliminated 
by brute force. Usually, for physical reasons, their function in statistics is successfully 
mimicked by their non-analytic Calogero substitutes ([]). 

We may summarize that the reality of spectra of our new, non-Hermitian and 
complexified Calogero model can be again tentatively interpreted as stemming from 
its multiple permutation VT symmetry. We can conjecture, therefore, that this type 
of a complex permutation could presumably play the role of a certain "weakened 
Hermiticity" in the many-body physics. This generalizes the parity x time-reversal 
symmetry which is already comparatively well understood within the VT symmetric 
quantum mechanics of isolated particles [FJ. 
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Figure captions 

Figure 1. PT symmetric energies for A = 2 

Figure 2. Toy spectrum 

Figure 3. VT symmetric energies for A = 3 
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Figure 2 . Toy spectrum 




